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ABSTRACT 

We present angular momentum transport (torque) measurements in two recent experimental stud- 
ies of the turbulent flow between independently rotating cylinders. In addition to these studies, 
we reanalyze prior torque measurements to expand the range of control parameters for the exper- 
imental Taylor-Couette flows. We find that the torque may be described as a product of functions 
that depend only on the Reynolds number, which describes the turbulent driving intensity, and the 
rotation number, which characterizes the effects of global rotation. For a given Reynolds number, 
the global angular momentum transport for Keplerian-like flow profiles is approximately 14% of 
the maximum achievable transport rate. We estimate that this level of transport would produce an 
accretion rate of M/Mq ~ 10"^ in astrophysical disks. We argue that this level of transport from 
hydrodynamics alone could be significant. 

Key words, accretion - accretion disks - turbulence - methods: laboratory - hydrodynamics 

1. Introduction 

Astrophysical disks are ubiquitous in the universe, orbiting around a wide panel of massive objects, 
from black holes to stars. An essential ingredient of their dynamics is the radial transport of angu- 
lar momentum that governs the rate of material falUng onto the central object. Reahstic models of 
this process are impeded by the complexity of the disk dynamics, including a wide variety of pro- 
cesses such as general relativistic effects, self-gravitation, radiation, dynamics of plasma magnetic 
fields and turbulence. In particular, the contribution from hydrodynamics to the angular momentum 
transport in such astrophysical flows is not yet understood. 

A natural first-step in the construction of a robust model is to consider the minimal list of 
ingredients that can capture and reproduce the observed disk properties. Such a model has been 
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proposed several decades ago by |Shakura & Sunyaev] ( |1973[ ) and |Pringle| ( |1981[ ), and was given 
the name a disks. In that model, the radial transport of angular momentum is parameterized by a 
turbulent viscosity coefficient, that needs to be prescribed, or measured. In the past, there have been 
a few attempts to measure this coefficient in numerical simulations ( Dubrulle|1992| l. 



An alternative promising option, as recognized by Zeldovich ( 1981| l and Richard & Zahn 
1999| l, is offered by focusing instead on laboratory experiments. Indeed, it may be shown ( Hersant 



et al. 112005 1, that under simple, but well founded approximations, the equations governing an 



disk are similar to the equation of motion of an incompressible rotating shear flow, with pene- 
trable boundary conditions and cyhndrical geometry. This kind of flow can be achieved in the 
Taylor-Couette flow (see Fig.[T]l, a fluid layer sheared between two independently rotating, coaxial 
cyhnders. 

This remark motivated several experimental studies of the transport properties in the Taylor- 
Couette flow, with contradicting conclusions. Depending upon the community, the strength of the 
turbulent transport is quantified in various ways. The most basic quantification ( |Wendt|1933| l is the 
one through the torque T which is necessary to keep the inner cylinder rotating at a given velocity. 
The dimensionless version thereof is G - T/ {Lmidpv^), where p is the fluid density, v its kinematic 
viscosity and Lmid the effective length of the cylinders ( Paoletti & Lathrop|201 1 1. (Eckhardt et al. 
|2007b[ ) and ( |van Gils et al.|[2011b| l used Nuoj - G/Gi^m, where Giam is the laminar torque, to 
quantify the angular velocity flux from the inner to the outer cylinder, in order to highlight the 
analogy between Taylor-Couette and Rayleigh-Benard flow ( DubruUe & Hersant|2002[ Eckhardt 
et al.|2007b|a l. The engineering community prefers the drag-coefficient c/ = G/Re^, where Re - 



-^\r]Re2 - ^eil and ?/ - a/b is the radius ratio and Re 1/2 the Reynolds number of the inner/outer 
cyhnder. In the astrophysics community, the same information is often expressed in terms of the y6 
coefficient, defined by y6 = IGif - rjfRe^) - 2jfcfl{n{\ - rjf ), which can be interpreted as 
a dimensionless turbulent viscosity. 



By reexamining previous data obtained by Wendt (Wendt 1933 1 and Taylor (Taylor 19361, 
Richard & Zahn ( 1999 1 showed that the jS-parameter is on the order of Prz ~ 10~^ for Keplerian 



disks. A more complete analysis, including results of new experiments on a classical Taylor- 
Couette flow from |Richard| ( |200"T] l, led to the same conclusion ( |Dubrulle et al.|2005] l. In contrast, 
an original experiment, built in Princeton ( Ji et al.|2006 1 so as to minimize end effects, measured a 
jS coefficient through radial Reynolds stresses and obtained a value smaller by more than an order 
of magnitude /3p ~ 7 x 10"^. Furthermore, the Princeton experiments could not distinguish the 
angular momentum transport in Keplerian flows from their measurements in solid-body rotation, 
where the cylinders rotate together thereby resulting in zero angular momentum transport. 

In the present contribution, we report precise angular momentum transport (torque) measure- 
ments in two independent experiments performed at very large Reynolds number and compute the 
resulting /? coefficient for different rotation profiles. To further study the scahng of angular mo- 



mentum transport in rotating shear flows, we compare our measurements to those of Wendt (19331, 
|Taylor| ( [T936| and [Richard] ( |200T] l; |Dubrulle et al.| ( [2005l ). We find a universal scaling of the torque 
for various rotation profiles that captures the effects of the various geometries used in the experi- 
ments summarized here. 
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Fig. 1. Taylor-Couette flow may be used to study angular momentum transport in rotating shear 
flows. The fluid is sheared between two independently rotating cylinders of radii a < b that rotate 
at angular velocities Qi and Q2 for the inner and outer cylinders, respectively. Although the total 
fluid height is given by L, some experiments only measure the contribution to the torque over 
a central section of length Lmid to minimize end effects, as in ( Lathrop et al.||1992b|a[ Lewis & 
Swinney|[T999l [van den Berg et al.|[2003l [2005] [20071 [Paoletti & LaQirop[[20n] [van Gils et sT. 



|2011bD . 



2. Apparatus and experimental details 

2.1. Apparatus: generalities 

The experiments presented here all examine the transport of angular momentum (torque) in classi- 
cal Taylor-Couette flow. Figure[T]shows a schematic representation of our experiments. The fluid is 
contained between concentric cylinders of radii a < b. The inner (outer) cylinder rotates at an an- 
gular velocity of Qj (fla)- The overall height of the fluid contained between the axial boundaries is 
given by the length L. The geometry of a particular Taylor-Couette flow apparatus is often defined 
in dimensionless form by the radius ratio 77 = a/b and the aspect ratio T = L/(b - a). 

The ideal Couette flow is infinite in axial extent. As such, several methods have been employed 
to handle the finite-size effects that are present in any laboratory experiment. The most common 
case has the axial boundaries rotate with one of the two cylinders. In other experiments, the axial 
boundaries are divided at the mid-gap and have the inner (outer) portion rotate with the inner (outer) 
cyUnder. This split-ring approach has been further developed to allow for pairs of rings that can 
independently rotate with respect to the cylinders. Allowing the axial boundaries to independently 
rotate aids in the suppression of finite-size effects as detailed in ( Ji et al.[[2006[ Schartman et al. 
[2009) [Burin et al.|2010| l. 

The data that we analyze and present stems from several diff'erent experiments, which are sum- 
marized in Table [T] The early experiments by |Wendt[ ( |1933| l measured the torque required to drive 



3 



Paoletti et al. : Angular momentum transport and turbulence in Keplerian flows 

Table 1. A summary of the parameters used in the experiments discussed here. The geometry of 
the apparatus is given by the outer cyhnder radius b, the radius ratio t] - a/b and the aspect ratio 
r - L/{b - a). The range of turbulence intensities are governed by the inner and outer cylinder 
Reynolds numbers Rei and Re2- The effects of rotation are parameterized by the rotation number 
Rn introduced by |Dubrulle et al.| ( |2005] l. 



Experiment 


b (cm) 


'7 


r 


max(Rei) 


max(Re2) 


min(/?n) 


max(Rn) 


Wendt 


1933 


14.7 


0.68 


10.64 


8.47 X lO'* 


1.24 X 10^ 


-0.71 


0.47 


Wendt 


1933 


14.7 


0.85 


22.93 


4.95 X lO'* 


5.83 X lO'' 


-0.80 


0.18 


Wendt 


1933 


14.7 


0.935 


52.63 


2.35 X 10-* 


2.52 X lO'' 


-0.42 


0.11 


Taylor 


1936 


8.11 


0.79-0.973 


50-384 





1.08 X 10^ 


0.03 


0.27 


Richard 


2001 


1 5.00 


0.72 


27.18 


9.50 X 10-* 


1.30 X 10^ 


-1.43 


-0.94 


Richard 


2001 


1 5.00 


0.72 


27.18 


9.50 X 10* 


1.30 X 10' 


0.38 


0.58 


Ji et al. 


2006 


20.3 


0.35 


2.1 


1.20 X 10^ 


4.50 X 10' 


-1.67 


-1.08 


Maryland 


22.09 


0.725 


11.47 


3.70 X 10^ 


1.50 X 10* 


-2.04 


1.46 


Twente 


27.9 


0.716 


11.68 


2.00 X 10^ 


1.40 X 10* 


-0.74 


0.17 



the inner cylinder for three values of 77. The system had a free surface at the top, although we 
quote an effective aspect ratio defined using the axial length of the apparatus. The data from |TayIor| 
1936| l determined the scahng of the torque for several values of tj with the inner cylinder station- 



ary (Qi = 0). Richard's experiments ( |Richard|2001[[Dubrulle et al.|2005[ ) did not directly measure 
the torque. Instead, critical Reynolds numbers for various transitions were observed and the anal- 
ysis presented by |Dubrulle et al.| ( |2005| l was used to estimate the corresponding torque. We do not 
present data from the experiments of Ji et al.| ( [2006[ ), but rather use thek reported measurements as 
a basis for comparison, as detailed in Section |4] The two most recent experiments by Paoletti & 



Lathrop ( 2011 1 and van Gils et al. ( 201 lb I are further described in the forthcoming sections. 



2.2. Apparatus: the Maryland experiment 



The experiments presented by Paoletti & Lathrop ( 201 1\ , henceforth referred to as the Maryland 
experiment, are conducted in the apparatus constructed by Lathrop et al. ( 1992b|a i, which was 
modified to allow the outer cylinder to rotate independently. The outer cylinder is an anodized 
aluminum cylinder with the same physical dimensions as in Lathrop et al. ( 1992b|a| l; Lewis & 



Swinney ( 1999[ l, specifically b - 22.085 cm and a working fluid height L = 69.50 cm. The inner 
cyhnder is stainless steel with a radius a - 16.000 cm yielding a radius ratio rj - 0.7245 and an 
aspect ratio F - 11.47. The inner cylinder is rotated up to Qi/27r = 20 s ' while the outer cylinder 
may be rotated in either direction up to \Q.2/27r\ - 10 s"^ Both angular velocities are measured by 
shaft encoders and controlled to within 0.2% of the set value. 

The axial boundaries in the Maryland experiment rotate with the outer cylinder To reduce end 
effects in the torque measurements the inner cylinder is divided axially into three sections of length 
15.69, 40.64 and 15.69 cm (see schematics in Lathrop et al. ( 1992b|a 1). Only the central section 
of length Lmid = 40.64 cm (see Fig. [TJ senses the torque of the fluid as described in |Lathrop| 
et al. ( 1992b|a i. The regions 2.58 gap widths from each of the axial boundaries, where secondary 
circulation setup by finite boundaries is strongest, are avoided in the torque measurements. 
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The desired accuracy of our measurements requires that the temperature of the water be pre- 
cisely controlled. In contrast to prior experiments ( Lathrop et al.|1992b|a Lewis & Swinney|1999 



[van den Berg et al.|[2003l [2005] [20071 |van Gils et al.|201 lb| ) where the system was cooled at the 
axial boundaries, we control the temperature of the water by heating and cooling the outer cylin- 
der. This allows the working fluid to be temperature-controlled along the entire axial length of the 
experiment, yielding a 6.5 fold increase in the temperature-controlled surface area. This is particu- 
larly important for Rayleigh-stable flows, where mixing is greatly reduced. In our experiments the 
temperature is 50 + 0.02°C yielding a kinematic fluid viscosity of v = 5.5 X 10"^ cm^/s, except for 
Re>2x 10*' where T = 90 °C and y = 3.2 x 10"^ cm^/s. 

2.3. Apparatus: the Twente experiment 

The Twente turbulent Taylor-Couette (called T^C) facility, here referred to as the Twente experi- 



ment, is described in great detail in |van Gils et al. (201 la i. In short, the working fluid height is 



L = 92.7 cm, has an inner radius of a = 20.0 cm and an outer radius of b - 27.94 cm. The maxi- 
mum inner and outer angular velocities are Q.\l2n - 20 Hz and |Q2/27!-| = 10 Hz, respectively. The 
system is fully temperature controlled through cooling of the upper and lower end plates, which 
co-rotate with the outer cylinder The system is operated with water at a temperature of 20 °C, 
resulting in a kinematic viscosity of v = 1 .04 x 10"^ cm^/s. The torque is measured over the middle 
part of the inner cylinder of height Lmid = 53.6 cm to minimize the influence of the end-plates 
(similar to ( Lathrop et al.|[T9"92b|a[ l). The torque is measured by a load cell imbedded inside the 



middle section of the inner cylinder in both the Twente and Maryland experiments, in contrast to 
measuring the torque on the drive shaft that would be affected by mechanical seals in the end plates 
or by velocimetry measurements, as in the Princeton experiments ( Ji et al.|2006 Schartman et al. 
[2009)[Bmm et al.|2010||Schartman et al.|201 l| l. 

2.4. Control parameters 



The analysis of our results is simplified by a proper choice of the control parameters. As in DubruUe 



et al. ( 2005 i, we use the foUowing three dimensionless control parameters, which fully describe the 



geometry and specify the state of the dynamical Taylor-Couette system: 

Re = -^\T^Re2-Re,\, (1) 
1+77 

^„ . (1_,)J?£1±J^, (2) 
riRe2 - Re I 



Rc = (3) 



where Rei = aQ.\dlv and Re2 - bQ.2d/v are the Reynolds number of the inner and outer cylinder 
and d - (b - a)is the gap width. 

The above control parameters have been introduced so that their definitions apply to rotating 
shear flows in general and not strictly only to the Taylor-Couette geometry. It is very easy in this 
formulation to relate the Taylor-Couette flow to the shearing sheet (plane Couette flow with rota- 
tion), by considering the limit Rq — > 0. The linear stability properties of the fluid can also be neatly 
recast using these control parameters. In the inviscid limit (Re — > 00), and for axisymmetric dis- 
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Table 2. A summary of the Taylor-Couette control parameters related to some astrophysical disks, 
as computed by 



Hersant et al. 



(2005 I. Here, b is the disk's outer radius, and M the accretion rate. 



Central object 


b (a.u.) 


1 


r 


max(Rei) 


rmx(Re2) 


Rn 


M/Mo 


TTauri 


1000 





0.94 


1012 


10'^ 


-1.33 


10-5 - 10-2 


FuOri 


5 





0.94 


10" 


10'^ 


-1.33 


0.1-1 



turbances, the linear stability properties of the flow are governed by the Rayleigh criterion, namely 
that the fluid is stable if the Rayleigh discriminant is everywhere positive: 

-d,L(r) > 0, (4) 
r 

where L(r) = r^Q(r) is the specific angular momentum. Applying this criterion to the laminar 
profile leads to 

(Rn + l)(Rn + 1 - ab/r^) > 0. (5) 

Since ab/r^ varies between l/rj and rj, one obtains that in the inviscid limit, the flow is unstable 
against infinitesimal axisymmetric disturbances when R'^^ < Rn < Rq^^ where = -1, re- 
spectively R"^^ - l/?7 - 1, are the marginal stability thresholds in the inviscid limit (superscript 
oo) in the cyclonic case (Rn > 0, subscript +), and anticyclonic case (Rq < 0, subscript -). This 
means that, in the anticyclonic case, any flow such that Rq < -1 is linearly stable against axisym- 
metric disturbances. However, it may be non-linearly stable or unstable against non-axisymmetric 
disturbances, as found in the experiments described here. 

In the present paper, we determine the intensity of the momentum transport by measuring the 
torque required to drive the inner cylinder for a given set of parameters {Re,Rn). Measurements 
described here suggest that at large enough Reynolds number, the nondimensional torque G - 
T /(pv^Ljnid), where T is the torque, p the fluid density, v the kinematic viscosity and L,„ici is defined 
in Fig. [T] varies approximately as G ~ Re" with 1.75 < a(Re) < 1.85, even for Rayleigh-stable 
flows. While this scaling exceeds the laminar case (G ~ Re), we note that these exponents are 
less than a simple dimensional argument would give: Such an argument would say that in a fully 
turbulent regime the transport properties could not depend on the molecular viscosity, which would 
imply that the quantity Cf ~ G/Re^ is Reynolds number independent. The deviations from G ~ Re^ 
mean that the viscosity is always relevant for the transport, due to its dominance in the boundary 
layers. Only when the boundary layers are sufficiently disturbed by introducing wall roughness 
does one obtain G ~ Re^, as we have shown in our previous work (van den Berg et al. 2003 i. 



2.5. Astrophysically relevant quantity 

In the astrophysical context, one often considers asymptotic angular velocity profiles of the form 
Q(r) ~ r^'' where q characterizes the flow. In that case q - -dlnQ/dln r - -2/Rn, which relates 
astrophysical profiles in the control parameters space of Taylor-Couette flows. For Keplerian flow, 
q - 3/2 resulting in Rq = -4/3, which is an example of a linearly stable, anticyclonic flow. While 
Keplerian profiles cannot be precisely achieved by Taylor-Couette flows, states with Rn < -1 are 
all linearly stable, anticyclonic flows that may be used to approximate astrophysical settings. We 
will refer to cases with 7?n < -1 as quasi-Keplerian, as in Ref. pi et al.|2006[ ). 
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Astrophysical disks are very wide with a <^ b. Therefore, their geometrical global value of 
Rq — > oo. In a local description of disks, one often considers the shearing sheet approximation, in 
which ?7 — > 1, so that Rq — > 0, as in plane-Couette flow. Here, we focus on global properties and 
consider only the limit rj ^ Q applicable to disks in the discussion. As a result of their astronomical 
size, the Reynolds numbers associated with accretion disks are very large. To give a few specific 
examples, parameters for typical TTauri and FuOri stars have been computed by Hersant et al.| 
( 2005 I and are given in Table [2] 

In disks, angular momentum transport intensity is quantified through the accretion rate M. In 
fact, as was shown in ( DubruUe et al.|20()5 Hersant et al.|2005 1, this rate is directly related to the 
turbulent viscosity and the dimensionless torque. This concept may be understood in the following 
way: pressure gradients substantially balance centrifugal forces that develop in (incompressible) 
laboratory flows. Keplerian disks, however, balance the centrifugal force with gravitational attrac- 
tion. Thus when angular momentum is transported radially outward, accretion disk matter responds 
by falling radially inward at the rate M. The flux of angular momentum, which is proportional to 
G, is thereby related to the accretion rate M by ( Dubrulle et al.|2005 Hersant et al. 2005| l, 
G n M 

where y6 is a turbulent viscosity parameter such that v, = fiSp-, with S a velocity shear, and Mq is 
an eff'ective accretion rate given by the mean surface density £, the disk inner radius a, mean radius 
? - {a + /7)/2,the mean angular velocity Q. and the mean disk height H through 

Mo = ta?Q.{^\ . (!) 



As such, measurements of G/Re^ in laboratory experiments may be used to estimate the hydrody- 
namical contribution to angular momentum transport in accretion disks. 



3. Results 

We study the scaling of the dimensionless torque G over a range of Reynolds numbers for various 
constant rotation numbers (see Fig. |2]l, spanning a wide variety of flow states, including linearly 
stable cyclonic and anticy clonic flows as well as linearly unstable states. Corresponding torque 
measurements from the Twente experiment are shown as red squares in Fig. [3] while the data from 
Maryland are represented by black circles. 



3.1. Existence of turbulence 

In laminar flow, angular momentum is transported only by molecular viscosity, so that the dimen- 
sionless torque G is proportional to Re, with the prefactor depending upon Rq. In Fig. [3] we show 
the expected laminar angular momentum transport (blue line) for the range of Re explored in this 
study for the case Q2 - 0. As can be seen all the torque measurements considered differ from the 
theoretical expectation for laminar flow in two ways: (/) they increase faster than Re and (//) they 
are much higher than the theoretical laminar value. This is an indication that the fluid is turbulent, 
with enhanced angular momentum transport in all the cases reported here. Note that these cases en- 
compass several situations with < - 1 (hnearly stable anticyclonic flow) and Rq > 0.4 (linearly 
stable cyclonic flow) in addition to many unstable states. 
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Fig. 2. Our experiments span a large range of the Reynolds number Re, which controls the turbu- 
lence intensity, and the rotation number Rq,, introduced by DubruUe et al. ( Dubrulle et al.|2005 l, 
which compares shear to overall rotation. The rotation number can be used to distinguish be- 
tween Rayleigh-stable flows {R^ < -lor/?n < (1- t])/t]) and those that are linearly unstable 
(-1 < Rn < {\ — rf)lTj). The data are taken from separate experiments by van Gils et al. ( 201 Vo) 
(red squares) and ( Paoletti & Lathrop|201 1| | (black circles). 



In particular, (approximate) Keplerian flows at Rn - -4/3 are turbulent at Re ~ 10^ in our 
experiments. This is in agreement with direct, visual observations of Richard ( 2001[ l, performed 
at somewhat lower Reynolds number. We argue that this indicates that quasi-Keplerian flows can 
efliciently transport angular momentum in spite of their linear stability. However, the nature of this 
nonlinear instability remains unclear. We suggest that systematically perturbing quasi-Keplerian 



flow states at lower values of Re while measuring the torque, as in the experiments by Taylor 
1936| l, could aid in the description of this likely nonlinear instability. 



3.2. Enhancement of momentum transport 

We now try to quantify the enhancement of angular momentum transport with respect to its laminar 
value, as a function of both Re and Rq. As in |Paoletti & Lathrop| ( |201 l[ ); [van Gils et"aL] ( |201 lb) , the 
torque at a given Re may be either increased or reduced depending upon Rn', the Reynolds number 
alone is insufficient to describe the transport. 



3.2.1 . Variation witli Rn 



Paoletti & Lathrop (201 1 1 observed that variation of the torque at a given Re shows a pronounced 



maximum as a function of Rn. An analogous dependence of the amplitude of an effective power- 
law scaling was also observed by van Gils et al. ( van Gils et al.|2011b l. To determine the effects 
of global rotation, we normalize G for each Re by G(Re, Rn - 0) ( [Dubrulle et al.|[2005| , which 
we denote as Gq. The dependence of G/Gq on Rn is shown in Fig.[5]with the data from Maryland 
shown as black circles and Twente in red squares. 

The agreement between the Twente and Maryland measurements is an important check, but 
not surprising given that the experiments have very similar geometries, measurement techniques 
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10' 



—I 1 1 1 — I — I — I — r- 



—I 1 r- 




Re{W) 



Fig. 3. The scaling of the dimensionless torque G with the Reynolds number Re is independent of 
the rotation number however the ampUtude varies as in Refs. ( [Paoletti & Lathrop|2011 [[van Gils| 
et al.|20lTb i. The measurements by Paoletti & Lathrop ( 201 l| i are in black circles while the data 
from [van Gils et aL[ ( [201 lb| l are shown by red squares. The thick blue line indicates the theoretical 
value for purely laminar flow, where angular momentum is transported only by molecular viscosity. 



and control parameters (see Table [T]). Therefore, it remains unclear if the scaling of the normalized 
torque with is universal or if it depends upon other parameters, such as the radius ratio 77. While 
the Twente experiment is capable of examining other values of 77 in future measurements, we may 
compare our results to previous experiments by Wendt ( 1933 1, Taylor ( 1936| l and Richard (2001 1, 
which were analyzed by [Dubrulle et al. (2005 1. The parameters of these past experiments are also 
summarized in Table [T| The measurements of Wendt (triangles), Taylor (orange diamonds) and 
Richard (magenta squares) are also shown in Fig.|4|a). 

As can be seen in Fig. |4] the dependence of the normalized torque G/Go on /?n for various 
radius ratios in the range 0.68 < rj < 0.973 collapse well. The data appear to follow distinct 
approximate scaUngs given by 



f(Rn) 



0.18 + 0.06 : Rn<-i 

1.79/?n + 0.05e-'*-^^'^" + 1.59 : - I < Rn < -O.IO 

-7.357?n + l-49e^-2'**" -0.51 : - 0.10 < /?n < 0.32 

-0.13/?n + 0.26 :7;n>0.32 



(8) 



Given the rather different span of Re in each experiment and the varying measurement techniques, it 
is not surprising that a few outliers remain. For example, the values of G/Go from Richard's exper- 
iments ( [Richard [[200 l| l (magenta squares) were computed using theoretical predictions regarding 



the correlation between the torque and critical numbers for stability, as described in DubruUe et al. 
( |2005l l. The data from |Wendt[ ( [T933] l, [Taylor] ( fT936| , [Paoretti & Lathrop| ( |20n] i and [van Gils etal^ 



(201 lb I are direct measurements of the torque, albeit with different apparati and measurement 
techniques. 
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1.5r 

1.25- 




Fig. 4. (a) The dimensionless torque G normalized by Go = G{Re,Ra - 0) for different values of 
the radius ratio tj nicely collapses when plotted as function of Rq. This normalization is different 
from Fig. 16 in Ref. (DubruUe et al. 2005 [ l where G is normahzed by Goo = G(Re, Q.2 - 0). 



The data from Wendt ( 1933 1 are shown as cyan (77 = 0.68), blue (77 = 0.85) and green (j] - 0.935) 



triangles. Additional radius ratios are examined using the data from Taylor ( 1936 1, shown as orange 
diamonds, Richard ("Richard 2001 Dubrulle et al. 2005| l, given by magenta squares, Paoletti & 
Lathrop ( 2011) indicated by black circles and van Gils et al. ( 201 lh\ by red circles. The solid lines 
are fits given by Eq. ([8| with smoothened transitions between the four different regimes. Keplerian 
flow states are indicated by the dashed line at Rq = -4/3. In (b) we focus on the peak of the curve 
shown in (a), but only present the data from Twente ( |van Gils et al.|[2011b I (red) and Maryland 
( [Paoletti & Lathrop|201 1\ (black). The solid line is the fit 



3.2.2. Variation witli Reynolds number 

The fits given in Eq. (j8]l may be used to compensate the torque to account for the effects of global 
rotation. Figure |5] shows the dimensionless torque G compensated by the fit f(Rn) given in Eq. (|8]l 
as a function of Re. The data from Twente are shown as red squares while those of Maryland are 
given by black circles. In both cases, the data agree well with the best fit given by |Lewis & Swinney] 
( 1999| l for the case Q.2 =0 (blue line). This indicates that the amplitude of the torque is affected by 
global rotation, as measured by G/Go, while the scaling with Re seems to be independent of Rn in 
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Fig. 5. The dimensionless torque G may be compensated by the fit f{Rn) given by Eq. The 
data from both Paoletti & Lathrop (201 1 1 (black circles) and van Gils et al. ( 2011b| l (red squares) 



agree with the best fit from Lewis & Swinney ( 1999 1 for the case where the outer cylinder was 
stationary (blue line). 



the regime of our investigation. Therefore the torque may be approximately factorized and written 
as 



G = f{Rn)8(Re), 



(9) 



where f(Rn) is measured in Fig.|4]and fit by Eq. ([8| and our measurements of g(Re), shown in Fig. 
|5] agree with previous measurements, and may be fit by effective power-laws or logarithmically 
varying functions such as the Prandtl-von Karman form (Lathrop et al. 1992b|a Lewis & Swinney 
[T999] l. 



3.2.3. Comments 

In addition to the observed data collapse, other interesting features arise in the normalized torque 
scaling with Rn- The data from Taylor's experiments ( Taylor|1936 1 were obtained with a stationary 
inner cylinder (Qi - 0), and the flow is therefore Rayleigh-stable. However, the data from |Weiidt| 
( [T933] l, Maryland ( [Paoletti & Lathrop|[20n) and Twente ( |van Gils et~aL]|201 lb i over the same 
range of the rotation number, namely < Rq < 0.27, are obtained with counter-rotating cylinders. 
This configuration is linearly unstable, yet the measured values of G/Go agree between all four 
experiments. We argue, as done by |Taylor| ( |1936| l, that the flow states in Taylor's experiments 
could not have been laminar and must have undergone a subcritical transition to turbulence. 

[Paoletti & Lathrop] ( [201 1[ ) observed that the torque scaling varied linearly with Rn in four 
regions of the parameter space corresponding to: region I: Rn < -1, region II: -1 < Rq < -0.10, 
region III: -0.10 < Ra < 0.38 and region IV: 0.38 < Rn, where we have reordered the region 
numbers to correspond to increasing Rq. Figure]?] shows that all of the data described here follow 
the region I-region II and region Il-region III crossovers. The crossover between regions I and II 
corresponds to the Rayleigh stability criterion for all values of 77. The crossover between regions II 
and III occurs at the observed maximum in the normalized torque G /Go as a function of Rq. While 
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Fig. 6. The scaling of (a) G/Re^, which is proportional to the accretion rate M and viscosity param- 
eter yS, decays slowly with Re while having an amplitude that varies with the rotation number Rq 
(see Fig.|4|l. (b) The effects of global rotation can be accounted for by compensating G/Re^ by the 
fits f(Rci) given in Eq. (|8]l. The measurements from Twente are shown as red squares. The Maryland 
data are given by black circles, except for quasi-Keplerian flows (Rq < -1), which are indicated 



by blue diamonds. The green continuous line is the theoretical fit given in Eq. ( 10 1 normahzed by 
fiRn = -0.2755), which corresponds to Q.2 - 0. 



this maximum occurs for all of the Wendt ( T933JI, Maryland (I Paoletti & Lathrop|201 l| l and Twente 
( van Gils et al.|2011b i datasets, there is no general theoretical prediction for the exact location of 
this maximum yet. 

The last crossover observed by Paoletti and Lathrop, between regions III and IV, corresponded 
to Qi =0 with Qi > and Q.2 < in region III (counter-rotation) to f2i > and Q.2 > in 
region IV (co-rotation). It is unclear whether this crossover is set by Rn = 0.38 or whether it 
always occurs at /?n(Qi - 0). The data from Taylor ( 1936[ ) has Qi = but given the values of t] 
used in those experiments Rq < 0.38. The Richard data for Rq > 0.38 fall below the observations 
by Maryland, although the scaling with Rq is similar It is unclear if this disparity is a result of the 
analysis used to deduce the torques from Richard's measurements or the difference in the control 
parameters from the Maryland experiment (see Table [T]). It would be fruitful to have independent, 
direct measurements of the torque in region IV to compare against those of Maryland at different 
radius ratios to better understand this final crossover. 

The Richard and Maryland data both indicate that the observed torque for quasi-Keplerian flows 
states (Rn < -1) is approximately 14% of the maximum observed torque for a given Reynolds 
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number. This is in contrast to the results of |Ji et ar] ( |2006[ l, who deduced that such flow states 
are "essentially steady" and unable to transport angular momentum hydrodynamically. In addition 
to the diff'erent geometries between the experiments in |Ji et al.| ( |2006| l and those of Richard and 
Maryland (see Table [TJ, the measurement techniques and finite-size eff'ects were different between 
the experiments. This contrast is further discussed in Section |4] Clearly, additional independent 
measurements of the angular momentum flux for quasi-Keplerian flows could aid this debate, as 
recently discussed by Balbus ( 2011| l. 



3.2.4. Astrophysical implications 



As described in Section 2.5 the scaling of the torque with Reynolds number may be used to de- 
termine either the turbulent viscosity parameter y6 or the accretion rate M through Eq. Our 
measurements of G/Re^, which is proportional to both j3 and M, are shown in Fig.[6[a). The mea- 
surements from Twente are shown as red squares while those of Maryland are given by black 
circles. The data indicated by blue diamonds are quasi-Keplerian flow states (Rq < -1) from the 
Maryland experiment. The quantity G/Re^ decreases slowly with Re (only logarithmically, see 
( DubruUe et al.|2005 Hersant et al.|2005 Grossmann & Lohse|201 l| i) as expected since G oc Re" 
with 1.8 < a{Re) < 2.0 for 10^ < Re < oo ( jLaflirop et aL]|l992b|a[ [Lewis & Swinney||l999| ). 
The amplitude of G/Re^ for a given Re is determined by Rq. The eff'ects of global rotation may be 
accounted for by compensating the measurements of G/Re^ by /(Rn) (see. Eq. (jsjl), as shown in 
Fig.|6|b). As discussed in Dubrulle & Hersant ( 2002[ l; Hersant et al. (2005 1, the scaling may be fit 
by a logarithmically decaying function, such as: 



1 



Re^ " (1 - 7j)y^ \n[j]Hl - rfjRe^lK^VI^-' 
with the values chosen to be K-j - 0.4664 and /Tg 



(10) 



10 to describe the best fit for the torque scaling 



provided by Lewis & Swinney ( 1999 1 for the case Q2 = 0. An alternative logarithmic dependence 
between G/Re^ and Re has been suggested by Grossmann & Lohse (201 1 1 (for the analogous case 



of ultimate Rayleigh-Benard flow), and it equally well fits the data. 

The expected accretion rate, measured by G/Re^ = 3.3 x 10"^, for quasi-Keplerian flows (blue 
diamonds in Fig. [6]) is approximately 14% of the maximum accretion rate, which occurs for Rq - 
-0.10. This corresponds to a value of /SmcI = 2 x 10"^ or M/Mq = 3.3 x 10"^^ at /?e = 7 x 
10^. Extrapolating towards astrophysical Reynolds numbers by using the formula of Hersant et al. 
( 2005| l, we get a value of M/Mq = 10"^ that is consist with values observed in disks around TTauri 
stars. 



4. Discussion 

The measured accretion rate for all of our quasi-Keplerian flows (Rq < -1) is approximately 14% 
of the maximum accretion rate, which occurs at Ra - -0.10, resulting in j0Md = 2 X 10"^. This is 



close to the value computed by Richard and Zahn /3rz = 10 ^ and to the value of Dubrulle et al. 
(2005 I ySo = 8 X 10"^. This is somewhat surprising given that quasi-Keplerian flow profiles are 



linearly stable while the maximum accretion rate occurs when the cylinders counter-rotate, which 
produces destabilizing shear that dominates the dynamics. This result is in stark contrast to the 
results presented by | Ji et aL] ( |2006| l and |Schartman et aL] ( |201 l| l, where /3p was measured to be an 
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Fig. 7. Finite-size effects may be characterized by axial profiles of the angular velocity Q. measured 
at the mid-gap (r - (a + b)/2) in the Twente experiment. In this example, the cylinder Reynolds 
numbers are Rei = 1.0 x 10^ and Re2 = (outer cylinder stationary). The global torque measure- 
ments are taken for 0.22 < z/L < 0.78 with the top boundary of these measurements indicated 
by the dashed line at z/L = 0.78. The resolution of the laser-doppler anemometry prohibits mea- 
surements of the angular velocity in the boundary layer near z/L = 1 where Q. - Q.2 - 0, m this 
case. 



order of magnitude smaller, specifically jSp = 7 x 10 ^. Thus, they concluded that hydrodynamics 
alone cannot efficiently transport angular momentum in quasi-Keplerian flows. 



Schartman et al. (2011 1 state that these incompatible results are caused by differing interpre- 



tations of finite-size effects, namely Ekman circulation produced by the axial boundaries. In the 
Princeton experiments ( Ji et al.|[2006 Schartman et al.|2009 Burin et aL]|2010 Schartman et al. 



201 l| l the aspect ratio is small (F x 2) and therefore finite-size effects would likely dominate. To 



minimize Ekman circulation, their axial boundaries are split into pairs of rings that can indepen- 
dently rotate with respect to the cylinders. This method has been shown to reduce the effects of 
Ekman circulation in the bulk of the fluid ( [Schartman et al.|2009| l. Furthermore, when Ekman cir- 
culation is not minimized their measured values of /3 for quasi-Keplerian flows are approximately 



half of the typical values presented here. Thus, Schartman et al. ( 201 1 1 argue that our larger values 
of yS are the result of Ekman circulation that is present in other experiments, but is minimized in 
their apparatus. 

We claim, on the other hand, that our measured torques cannot be solely attributed to finite-size 
effects. The main premise of our argument is that the normalized torques G/Gq agree quite well for 
all of the different experimental data detailed here, even though the apparatus geometry, measure- 
ment techniques and control parameters vary greatly. Given this widespread disparity in geometry, 
rotation rates and range of Reynolds numbers, one would expect that the effects of Ekman cir- 
culation would differ and result in systematic discrepancies between the experiments. We do not 
observe such discrepancies and we further detail our arguments along these lines below. 

The torques in the Maryland and Twente experiments are measured over only 56% of the axial 
length of the flow centered at the mid-height of the apparatus. This design has been intentionally im- 
plemented such that the torque measurements are unaffected by flows within 2.58 radial gap widths 
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of either axial boundary where Ekman circulation is strongest. Wendt's experiments ( |Wendt|1933| l 
had a free upper surface, which is therefore devoid of Ekman circulation. Even though Taylor's 
measurements ( |Taylor|1936| l were affected by Ekman circulation for F > 100, his measurements 
of the torque for Rayleigh-stable flows that are expected to be most affected by Ekman circulation 
agree with our measurements of counter-rotating flows for the same value of Rq, where Ekman 
circulation is dwarfed by the dominant shear 

One of course wonders on whether this Ekman circulation caused by the co-rotating upper and 
lower plates affect the flow patterns. For laminar flow e.g. for purely outer cylinder rotation this 
clearly is the case ( [Coles & van Atta|1966[ |Hollerbach & Fournier|2004| . However, we argue that 
the strong turbulence destroys such an axial dependence of the flow field. To check this assumption, 
we have measured the full turbulent velocity profile with Laser Doppler Anemometry (LDA) ( |van| 
Gils et al.|2011c I. The result is shown in Fig. |7] Evidently, over the length of the central section 



of the inner cylinder where we measure the torque there is no visible z-dependence of the mean 
velocity; the plate boundary layers only produce an axial dependence far above the domain of the 
central section of the inner cylinder. 



Ji et al. ( 2006 1 stated that their measured velocity fluctuations and values of yS were much 



higher for more viscous fluids, even though the flow profiles were the same. They interpreted this 
observation to imply that the residual finite-size effects penetrate deeper into the bulk of the fluid 
at lower Re, even though Ekman circulation was claimed to have been minimized in both sets of 
measurements thereby resulting in similar (nearly ideal-Couette) flow profiles. The values of Re 
for the experiments presented here span two orders of magnitude while the values in viscosity vary 
by a factor of 34, yet no such systematic effects on the measured torques have been observed. 
Furthermore, given the small aspect ratio of the Princeton experiments one would expect Ekman 
circulation to produce significantly higher values of /3 when the axial boundaries co-rotate with 
the outer cylinder, as in the experiments presented here. However, even in this case their measured 
values are a factor of two smaller than those presented here and in Paoletti & Lathrop ( 2011| l. 

The normalized torque G/Go is approximately constant for each value of Rn characterized by 
the present experiments. This indicates that the ampUtude of the torque varies with Rq, however 
the scaling with Re is unchanged (see Fig. [5]l. [Paoletti & Lathrop | ( |2011| l observed that their torque 
scaling for Rayleigh-stable and unstable flows agreed with the Prandtl-von Karman description 
of [Lathrop et al.[ ( [l992a| l; [Lewis & Swinney| ( |1999| l. The Twente data, initially presented by [van| 
Gils et al. (2011bi, was well-described by the theoretical predictions of Eckhardt et al. ( 2007b[ l 



for all rotation numbers explored. If Ekman circulation significantly contributed to the torque for 
Rayleigh-stable flows, as argued by Schartman et al. ( 201 l| l, then one would expect systematic 
disparities between these theoretical descriptions and the measured torques. As shown in Fig.|5] 
though, the torque scaling agrees well with the observations of Lewis & Swinney ( 1999 1 for the 
case Q2 = 0/or all values ofRa. 

The Princeton experiments ( Ji et al.|2006 Schartman et al.|2009 Burin et al.|2010 Schartman 



[et al.[[20TT[ ) measure local velocity correlations to determine the global angular momentum flux, 
whereas the Maryland and Twente experiments directly measure the global angular momentum 
flux using precision calibrated force-arm measurements. The distinction between the measured 



quantities has made direct comparisons of the experiments difficult. To bridge this gap. Burin et al. 



( [2010| l recently used their local velocity measurements to determine the corresponding torque G. 
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To compare to previous studies with Q2 = ( [Lathrop et al.|1992b|a[|Lewis & Swinney|1999| ) |Burin| 



et al. 



(20101 fit the scaling of c^'^^ to Xog^^fjyRe * c/), where cy = GjRe^, and determined a slope of 
1.40 ± 0.16. The authors claim that this value agrees with the measured value of 1.52 reported by 
Lathrop et al. ( 1992b|a l and 1.56 measured in Lewis & Swinney ( 1999 1. However, the appropriate 



Prandtl-von Karman scaling predicts the following relationship ( Lathrop et al.|1992b|a Lewis & 
Swinney|1999) 



= A?logio(/?eV^) + M, 



(11) 



in contrast to the scaling shown in Fig. 5 of Burin et al. ( 2010| l. Furthermore, the fit parameters 
and M are predicted to depend upon the radius ratio in the following way 



N = [{\-rf)\n\Q]l{riKylln), 



M : 



InTo 



In 



1 



1 +'7/yJ V2^ 



+ KyQ 



(12) 
(13) 



where a: = 0.40 is the von Karman constant and -5. The radius ratio 77 = 0.35 in the Princeton 
experiments, therefore the Prandtl-von Karman skin friction law predicts N - 5.16 and M - -3.15. 
These values do not agree with our analysis of the data in |Burin et al.| ( |201^ where N ~ 2 and M ^ 
-0.4. In fact, the measured values of G determined by the local velocity fluctuations are on average 
five times larger than those predicted by the Prandtl-von Karman skin friction law that has been 
verified using direct torque measurements (Lathrop et al. 1992b|a[ Lewis & Swinney 1999 Paoletti 



|& Lathrop |201 l| l. Therefore, we argue that the local velocity fluctuation measurements used in the 
Princeton experiments ( Ji et al.|[2006 Schartman et al.|2009 Burin et aTjpOlO Schartman et al. 



|201 l| l face severe experimental challenges. In fact, |Huisman et aL] ( |2011| l showed that the local 
angular momentum transpoit undergoes fluctuations that are two orders of magnitude larger than 
the mean value. However, once the local angular momentum transport is averaged long enough and 



over a large enough region, Huisman et al. (2011 1 achieved perfect agreement between the local 
and the global measurements. 



Clearly the contrasting conclusions of the Princeton experiments ( Ji et al. 2006 Schartman et al. 



|2009[ |Burin et al.|2010[ [Schartman et al.|[2011| l and those presented here warrant further experi- 
mental and theoretical investigation, as suggested by Balbus ( 2011[ l. It may be that the Princeton 
measurements of the local velocity fluctuations may not accurately determine the global angular 
momentum transport (( Burin et al.|2010} ), due to lack of convergence. On the other hand, the ef- 
fects of the axial boundaries on the global torque measurements presented here have not yet been 



characterized in detail. The Twente experiments (Huisman et al. 201 1 1 provide both local velocities 
as well as the global torque, which thereby presents future opportunities to clarify these distinc- 
tions. Vice versa, direct global torque measurements in a split axial-ring apparatus, such as the 
Princeton device, would also be scientifically useful. 



5. Conclusions 



The measurements that we present here indicate that the dimensionless torque G may be described 
as a product of functions that separately depend upon the Reynolds number Re and the rotation 
number Rq. The rotation number Rq introduced by DubruUe et al. ( 2005 | l collapses all of the 
experimental data presented here for various radius ratios, aspect ratios and Reynolds numbers 
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(see Table [T]l. The effects of global rotation, measured by G/Go, where Go = G(Re,Rci - 0), 
are shown to be well-described by the function f{Rn) given in Eq. ([8]) for the entire range of Rq 
explored in the present studies. 

The most astrophysically relevant result is the scaling of G/Re^, which is proportional to both 
the turbulent viscosity parameter /3 and the accretion rate M. We observe that the expected accretion 
rate for quasi-Keplerian flows (Rn < -1) is approximately 14% of the maximum attainable rate 
for a given Reynolds number We therefore get a dimensionless angular transport parameter y^Md - 
2 X 10"^, in agreement with Richard and Zahn but more than an order of magnitude larger than 
the Princeton experiments. Such angular momentum transport is able to produce dimensionless 
accretion rates on the order of M/Mq ~ 10"^, which are compatible with observations in disks 
around TTauri stars ( [Hersant et al.|2005] l. We have argued that our results are not easily attributed 
to finite-size effects, such as Ekman circulation. As such, we conclude that hydrodynamics may be 
able to transport angular momentum at astrophysical rates in spite of the linear stability provided by 
the radially-increasing angular momentum profile. This level of transport must result from turbulent 
means, thereby implying that quasi-Keplerian flows can be nonlinearly unstable to finite amplitude 
disturbances. 
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